Simple method for excitation of a Bose-Einstein condensate 
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I. INTRODUCTION 



Soon after the first spectacular realizations of the Bose- 
Einstein condensate (BEC) in cooled and trapped atomic 
gases [f]-pl , investigations of possible new effects involv- 
ing the condensate appeared. The BEC allows us to 
study several typical quantum mechanical phenomena on 
a macroscopic level - because the macroscopic sample of 
atoms is described by a single wavefunction. A stan- 
dard example is the splitting of the condensate into two 
spatially separated parts y followed by a superposition 
of the parts. The observation of the interference fringes 
[Q-Q is a manifestation of the quantum coherence be- 
tween two macroscopic parts of the condensate. By leak- 
ing the atoms from the condensate (typically downwards 
- due to gravity) one may prepare an "atom laser" M,p| . 

It has been also realized that collisions between spa- 
tially separated condensates may be used to create col- 
lective excitations in the condensate |^,|l2l. Assume a 
standard mean field single particle description of the gas 
of weakly interacting bosons in the limit of vanishing tem- 
perature (for reviews see 0-03). The time-dependent 
equation governing the state of the BEC is then the cel- 
ebrated Gross-Pitaevskii equation (CPE). Its solutions 
may describe either solitary waves or vortices as is typ- 
ical for the nonlinear equation. For characterization of 
soiitons the language of nonlinear optics Q is quite use- 
ful. One may consider then bright soiitons (bell shaped 
structures propagating without dispersion) , dark soiitons 
(with a node in the middle - an analog of the first ex- 
cited state in the non-interacting particles picture) or the 
intermediate grey soiitons. 

Not only collisions may be used to create excitations 
of the BEC. In fact several schemes have been proposed, 
some of them being successfully applied in experiments. 
It has been suggested that a resonant Raman excitation 
scheme may be utilized to excite vortex states [ISl . How- 
ever, the resonance is modified appreciably during the 
process of transferring the population from the ground 
state to the vortex state due to the nonlinearity of the 
GPE. An apparently more robust approach is the adia- 
batic scheme of Ml which takes full account of the non- 



linearity. It utilizes a controlled laser induced adiabatic 
transfer, populating solitonic or vortex solutions of GPE, 
depending on the details of the excitation. The adiabatic 
transfer uses internal atomic transitions combined with 
appropriate states of the condensate. A phase imprint- 
ing method, originally proposed in [1^], produces a phase 
shift between two parts of the condensate. This method 
in fact has been applied experimentally to create dark 
soiitons both in cigar shaped BEC |jl^ and in the spher- 
ically symmetric condensate [|9|. 

Another technique based on laser stirring of the con- 
densate allows the production of several different vortex 
states [|0|. Recently the method has been applied to 
create lattices containing over 100 vortices pi[| . 

The aim of this paper is to discuss in detail yet an- 
other method which, in our opinion, may serve to gen- 
erate collective excitations in a BEC. It allows for the 
creation of grey (or even dark) soiitons as well as vor- 
tices. The method which we propose resembles to a cer- 
tain extent the adiabatic passage scheme of [Q. In the 
latter, the transfer of population between two internal 
atomic states is accompanied by an appropriate change 
of the condensate wavefunction into a dark soliton, two- 
soliton or vortex solution of the GPE [|l6[ . Our method 
originally proposed for non-interacting particles 122] and 
extended to weakly attractive interaction assumes a fast 
sweep of the laser beam across the trap. In this way the 
trapping potential varies with time enabling a transfer 
of population to excited BEC states. Under appropriate 
conditions this approach allows for an efficient creation 
of collective excitations in the condensate as explained 
below. 

In Section II we present the method for a non- 
interacting particles model. We discuss both the exci- 
tation of soiitons in the one dimensional (ID) case |2^ ] 
as well as the possibility of vortex creation in the effec- 
tive two-dimensional (2D) example. In the next sections 
we extend the approach to interacting particles. In Sec- 
tion III we apply it to particles with repulsive atom-atom 
interactions (positive scattering length ao) in ID model 
case. In this way we complement our earlier study for 
attractive interactions p^ for the case most often met 
in experiments with a BEC |l|,^. Extending the treat- 



merit to 2D, we also consider the excitation of vortices 
for interacting particles. In Section IV we show some 
related observations on behavior of GPE energy levels 
while changing a parameter of the potential. 



II. NON-INTERACTING PARTICLES 

Let us first consider the simplest situation: The case 
of non-interacting particles. A BEG of non-interacting 
particles is not realized in nature, yet it may serve as 
a good model to describe the basic idea underlying the 
scheme we proposed. For an excitation of solitons prop- 
agating along a given direction ID model is clearly suf- 
ficient. For excitation of vortices, though, at least a 2D 
model is required. Such low-dimensionality models may 
be fully justified for non-interacting particles, assuming 
separability of the trapping potential. For interacting 
particles such simplified models may be of value for ap- 
propriately prepared (cigar shaped or flat disc shaped, 
respectively) condensates pl]- p^p3| -p7[ . 

Let us discuss the excitation of solitons first. Con- 
sider the condensate which occupies the ground state of 
the harmonic trap. Let us now sweep the region where 
the condensate is located with an additional laser beam, 
whose frequency is appropriately tuned close to the reso- 
nance of some internal atomic transition. By an adiabatic 
elimination of the upper atomic state (possible if the laser 
is detuned from the exact resonance) one may show p^ ] 
that such a laser beam creates an effective additional po- 
tential well (or a barrier, depending on the sign of the 
detuning with respect to the atomic transition) for the 
motion in the external atomic degree of freedom. The 
laser intensity is typically Gaussian-shaped in the direc- 
tion perpendicular to the direction of propagation. The 
atoms experience a potential proportional to the inten- 
sity which may be represented as 
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We use the trapping harmonic oscillator units, i.e. tot for 
time and ^fijmu! for length, where uj is harmonic oscil- 
lator frequency while m stands for atomic mass. (These 
units are also used through the rest of the paper). A sim- 
ilar modification of the potential has been used already 
in experiments to split the condensate into two parts [^ . 
We propose modifying the trapping potential in a dif- 
ferent way. Assume that the local Gaussian well is cre- 
ated on the very edge of the harmonic potential well thus 
not affecting the condensate. Then we change the laser 
beam direction slowly, in effect sweeping the well across 
the trapping potential. At the same time we gradually 
decrease the laser intensity, thereby decreasing the depth 
of the well. Such a procedure is equivalent to a change 
of xq from some negative value to zero in (Q) assuming 
Uq positive. 



The procedure proposed has a clear quantum mechan- 
ical meaning. For a sufficiently slow change of the poten- 
tial the system adiabatically follows its quantum energy 
levels except in the vicinity of avoided crossings. The 
energy gap of crossings may be controlled by choosing 
appropriate values of C/q and a in (^. In particular it 
is easy to arrange that the avoided crossing between the 
ground and the first excited state of the potential occurs 
during the sweep of the local potential well (see Fig. |^). 
Moreover such an avoided crossing may be made suffi- 
ciently narrow to be passed diabatically during the po- 
tential sweep. Then, when the local potential well disap- 
pears the particle, originally in the ground state of the 
harmonic trap, is left with a high probability, pi, in a 
first excited state (the Landau-Zener transition) . The ef- 
ficiency of the process depends on the size of the avoided 
crossing (which should be much smaller than the mean 
splitting between levels) and how quickly the potential is 
modified. 

As discussed in the earlier report [E2| , numerical simu- 
lations fully confirm the proposed scheme. Without any 
special optimization attempt, choosing the parameters 
of the potential (|^) as Uq = 13.4, a = 0.2 and chang- 
ing xq from —7 to with the velocity io = 0.02 yields 
Pi — 0.99. Similar values of pi are obtained for dif- 
ferent values of Uq and cr. The method is also robust 
with respect to the functional form of the potential well. 
We have checked that similar pi values are obtained if 
instead of the arctan(xo) we use other smooth monoton- 
ically changing functions of xq- 
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FIG. 1. Energy levels for a single particle in the potential 

m) for Uo = 13.4 and a = 0.2 as a function of xq- Note the 
narrow avoided crossing between the ground (the lowest solid 
line) and first excited (the lowest dashed line) states around 
xo ~ —4.5. Similar avoided crossings occur between the first 
and second excited states, between the second and third ones 
and so on. 

Looking at the level dynamics as a function of xq (com- 
pare Fig. 1^) one readily realizes that it is easy to general- 
ize the excitation mechanism to obtain higher excitations 



of the condensate. It is sufficient to sweep the local poten- 
tial well twice in order to get a highly efficient transfer of 
population from the ground state to the second excited 
state of the potential. This extension is easily tested 
numerically in the model with non- interacting particles. 
Keeping the same parameters as in the single excitation 
case but merely repeating the potential sweep the second 
time we get p2 = 0.99 as the squared overlap between the 
final wavefunction and the second excited state of the 
trap. A third consecutive sweep yields a "triply" exited 
state with the probability p^ = 0.99 again without any 
modification of the parameters of the sweeping potential. 
All these tests of single as well as multiple excitations 
indicate that, while we propose to excite the condensate 
by sweeping the trapping potential using the local po- 
tential well, the excitation may in fact be realized in a 
number of different ways. The key feature necessary for 
our method is the presence a narrow isolated avoided 
crossing between the ground and excited states. 



U{x,y) = 



x^ + y^ 




FIG. 2. Energy levels, in the rotating frame (SI = 0.6), for 
a single particle in the potential (pH, for Uo — 25 and a = 0.2, 
as a function of xo. Note the narrow avoided crossing be- 
tween the ground (the lowest solid line) and first excited (the 
dashed line) states around xo — —4.5. For xq = the en- 
ergy levels correspond to 2D harmonic oscillator states with 
1/2 = 0, Lz — 1, Lz — 2 and L^ — 3 from bottom to top. 

Our method to create vortices in a 2D model system 
is a natural extension of the former approach. It consists 
of diabatic transition through an isolated avoided cross- 
ing. Now we consider a BEG in a cylindrically symmetric 
2D harmonic trap. An additional laser beam produc- 
ing a local potential well is now rotating with frequency 
57 around the symmetry axis of the trapping potential. 
However, the energy is not conserved: Looking at the 
system in the frame rotating with a frequency f2 we may 
consider the eigenvalues of the Hamiltonian operator 
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+ U{x,y)-nL, 
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t/oVarctan(|xo|)exp ( — ^ ) . (3) 



For different distances |a;o| of the laser beam from the 
center of the trapping potential we can calculate energy 
levels of the Hamiltonian (0) . For appropriate values of 
Uq, a and f2, when changing xo from some negative value 
to zero, one can observe a narrow avoided crossing be- 
tween the ground and first excited states (see Fig. ||). 
The latter corresponds, for vanishing laser beam, to the 
first excited state of the harmonic potential with L^ = 1. 

For the same parameters as used in Fig. |[ we perform 
a time dependent numerical simulation. Starting from 
a ground state of the condensate and changing xq from 
—5 to with a velocity xq — 0.036, we obtain the first 
excited state with L^ = 1 with more than 99% accuracy. 
One may envision that after a single sweep we illuminate 
a condensate second time with a similar sweep. This, via 
a second avoided crossing (compare Fig. g) yields Lz = 2 
state with high efficiency (more than 99%). 

Provided, therefore, that the resulting picture is not 
modified strongly by the interaction of particles consti- 
tuting a BEG the proposed method should be able to pro- 
duce both solitons and vortices with quite high efficiency. 
The next sections describe the effect of the interactions 
on the proposed mechanism. 



III. COLLECTIVE EXCITATIONS IN BEC OF 
INTERACTING PARTICLES 



While the proposed method seems to be quite robust 
for non-interacting particles model its applicability to a 
BEG of interacting particles is far from clear. After all, 
the particle interaction necessarily changes the energy 
levels of the system. The proposed method relies on nar- 
row avoided crossings between levels when changing the 
parameter of the system - it is thus sensitive to the de- 
tails of level dynamics. It is not obvious whether the 
presence of the interaction between the particles will not 
destroy the proposed mechanism of the excitation. 

The partial preliminary answer has been given already 
in [g2[ , where we considered the effect of attractive atom- 
atom interactions present in a BEG of Li atoms [pi on 
the creation of solitons in a ID model. For such a con- 
densate the number of particles is not too big and the 
effect of atom-atom interactions on the behavior of the 
system is rather small. By a direct integration of the 
time-dependent GPE in ID 
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where 



(starting with the ground state of the condensate in the 
harmonic trap for g = —5) we were able to get a 97.5% 



population transfer into a collective state correspond- 
ing to the first excited state in the independent particle 
model. In (Q) g is a measure of nonlinearity and is pro- 
portional to number of particles N in the BEC. The value 
g = —5 taken for numerical simulation corresponds, for 
a trap used in [g|, to iV = 900 atoms in the condensate 
fraction - a typical number in Li condensate H . 



where H is given by Eq. (Q) . We get the time independent 
version by substituting ip{x,y,t) = exp{—i^t)ip{x,y) 
with /x being the chemical potential. The resulting time- 
independent equation 



H(p + g\ip\'^V = A*'/' 



(6) 
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FIG. 3. Single particle reduced probability densities of 

the condensate (i.e. solutions of the time-dependent GPE M) 
for g — 50) after the first [panel (a)] and second [panel (b)] 
potential sweeping. 

Clearly this is not the full story. The ID approach 
for interacting atoms is not exact, since nonlinearity cou- 
ples different degrees of freedom. Still, a one-dimensional 
approach based on the GPE is often used and may be 
justified for asymmetric traps [pl-fs 23-27 



More 



portantly most of the condensates realized in laborato- 
ries consist of particles with repulsive atom-atom inter- 
actions. Such condensates can easily hold about N = 10^ 
particles pd] p^ . Consequently the nonlinear term in (^) 
becomes much more important than for the attractive 
interaction. Thus the true test of our method requires a 
simulation for large positive g values. 

Fig. Sa shows the final wavefunction obtained for the 
g = 50 case by a numerical integration of Eq. (^ , taking 
as the initial state the ground state of the condensate for 
the same value of g. The parameters of the potential are 
Uo = 13.4 and a = 0.2, and the velocity of the sweep is 
xo = 0.6. The overlap of the wavefunction depicted in 
Fig. ga with the ideal "excited" state of the condensate 
is pi =0.98 at the end of the potential sweeping. 

Similarly successful is a double application of the po- 
tential sweep in order to obtain a "two-node" collective 
state of the condensate. We use the same parameters as 
above. In fact for a second sweep we start from the wave- 
function shown in Fig. ^a and make the second sweep 
identical to the first one. The results are depicted in 
Fig. gb. The overlap of the wavefunction at the end of the 
sweeping with the exact solution of the time- independent 
GPE is P2 = 0.82. 

Consider now the excitation of vortices in the 2D 
model. The time dependent GPE in rotating x, y co- 
ordinates reads 



^^ = Hi^ + g\^\^i^, 
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FIG. 4. Panel (a): single particle reduced probability den- 
sity of the condensate (i.e. solution of the time-dependent 
GPE dl) for g = 100, Uo = 25, cr = 0.2 and fi = 0.23) at 
the end of the potential sweeping where xo has been changed 
from —7 to with a velocity io = 0.35, see (H). Panel (b): 
the same as in the panel (a) but for g = 500, f2 = 0.12 and 
Xo going from —9 to with a velocity of 0.53. 

can be solved by the method of self-consistent field or 
by using the imaginary time propagation approach. The 
former is quite robust for the ID GPE but for 2D case it 
does not work efficiently. The latter is extremely effec- 
tive for the "ground state" of the condensate, however, 
its application to "excited" states becomes difficult (es- 
pecially for a trap without rotational symmetry) since for 
g ^ the stationary solutions are no longer orthogonal. 
The third possibility is to find the solution of Eq. (|^) by 
minimizing (0|0), where 



(5) 



= Hif^ 



-gW?^- 



■ flip, 



(7) 



with 



/2=((^ii/+5i(^n^). 



(8) 



while the changing parameter corresponds to a fictitious 
time |0|-|3|]. 



Indeed, starting with some initial function (p decomposed 
in a given basis, standard procedures of minimizing of 
multidimensional functions can lead to a desired solution 
provided the initial guess function (p is sufficiently close 
to the exact solution Lp. 

With solutions of the time-independent GPE at hand, 
we may prepare a BEG in its ground state, integrate the 
time-dependent GPE with the spiral-like potential sweep 
(given by (^ in the rotating frame), and compare the 
final wavefunction with the excited vortex-like solutions 
of the time-independent GPE. The resulting wavefunc- 
tions are presented in Fig. |^ for t; = 100 and g = 500. 
Their overlaps with the excited state corresponding to 
one quantum of the angular momentum (per particle) 
are 0.99 and 0.98 respectively. The non-ideal population 
transfer is responsible for a slightly asymmetric shape of 
the final wavefunctions. However, the method may be 
quite effective in generating vortex-like excitations in a 
BEG. Let us note that our approach resembles to some 
extent the stirring approach to vortex creation ||20(| . In 
that method the potential as a whole is rotated with a 
certain frequency fJ. In our approach a static, cylin- 
drically symmetric potential is supplemented, by an ad- 
ditional laser beam, with a narrow (with respect to the 
BEG dimension) structure moving along the spiral. Thus 
the physical picture of transferring the angular momen- 
tum to the condensate is a bit different in both cases. 
Naturally it differs also from the phase imprinting tech- 
nique [|7|. 



IV. BEHAVIOR OF LEVELS FOR INTERACTING 
PARTICLES 



The original physical picture of the excitation scheme, 
as discussed in Section II for the non-interacting parti- 
cles model, is based on the diabatic transitions between 
levels via narrow avoided crossings. The actual imple- 
mentation for interacting particles has been tested by 
numerical integration of the time-dependent GPE with- 
out resorting to the explicit changes of GPE levels with 
respect to a modification of the potential. To see whether 
the same picture may be invoked for the interacting par- 
ticles we have solved the ID time-independent GPE not 
only for the harmonic potential but also in the presence 
of the laser beam. Obtained energy levels Ig^l ^rs shown 
in Fig. ^ as a function of the position of the center of 
the laser beam. Observe the presence of loop-like struc- 
tures. Such a behavior is impossible in the case of a 
linear Schrodinger equation. Actually, it can be shown 
that the changes of levels' energies may be considered as 
a true Hamiltonian classical dynamics where the energies 
of levels play the role of positions of fictitious particles 
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FIG. 5. Energy (i.e. chemical potential |2£f|) levels of the 
condensate in the potential (nl) for g = 50, Uq = 13.4 and 
o = 0.2 versus a value of the parameter xq, i.e. the position 
of the laser beam. 
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FIG. 6. Solid lines: energy (i.e. chemical potential [E9|) 
levels of the condensate in the potential (|l|) for g = —1, 
Uo — 6.4 and a = 0.5 versus a value of the parameter a^o, 
i.e. the position of the laser beam. Specific points of the 
ground state level are indicated by letters, see text. Dashed 
lines are energy levels of the corresponding linear Schrodinger 
equation. 

To explain the appearance of the loops let us consider 
solutions of the time-independent GPE for weak attrac- 
tive particle interactions. Figure |6| shows the lowest en- 
ergy levels for both the non-interacting and interacting 
case with g = —1 [see (^]. 

In the non-interacting case avoided crossings result 
from changing the shape of the double well potential. To 
approximately predict energy levels in a double well one 
may consider each well separately. Then for an asymmet- 
ric potential, the resulting energies in each well are usu- 
ally considerably different. However, for certain shapes. 



the energies become degenerate. Taking the tunnehng 
between the wells into account, the level crossing changes 
into an avoided crossing (see Fig. ||). In the g — case, 
the change of the ground state energy (when going from 
the left to the right in Fig. O) corresponds to the transfer 
of probability density from the harmonic well to the local 
well created by a laser. At the point of the avoided cross- 
ing both wells are equally populated i.e. the ground state 
consists of equally weighted symmetric superposition of 
eigenstates of the right and left well. 




FIG. 7. Panel (a): plots of the potential (hi) (dashed line) 
and the effective potential (0) (solid line) corresponding to 
point A in Fig. M. Panel (c): probability density of the con- 
densate at the point A in Fig. g. Panel (b) and (d) are the 
same as in the corresponding previous panels but for point B 
in Fig. |. 

Analyzing the interacting particle case we may con- 
sider a solution of the GPE as if it is a solution of the 
Schrodinger equation with an effective potential 



Veff{x)^V{x)+g\^{x)\ 



(9) 



For g < 0, the interaction term deepens a potential well 
in which the probability density is localized (see Fig. 0). 
Consider the ground state of the harmonic trap and sup- 
pose that the laser beam approaches the center of this 
trap. At each position of the beam ccq, we can calculate 
approximate energies considering each well of the effec- 
tive potential (H) separately. The resulting energy in the 
left well is higher than the one in the right well when the 
probability density is situated in the right well (compare 
Fig. R). This holds up to the point indicated as A in 
Fig. y, where the lowest eigenenergy of left well treated 
alone approaches the energy of the condensate in the har- 
monic trap. At that point we can start populating the 
left well. However, transferring particles from the right 
to the left well breaks the energy balance between the 
two because of the density term in the effective potential 
(^ . The only way to restore balance is to lift the left well 
a little, i.e. to shift the center of the laser beam back to 
the left. 



Transfer of particles from the harmonic trap to the dip 
accompanied by shifting of the laser to the left can be 
carried on until nothing is left in the former. When this 
stage is reached (point B in Fig. ^) , all particles are local- 
ized in the left dip. Now this dip can be lowered without 
populating the right well. This will decrease the energy 
of the condensate which at some point matches the value 
of the starting energy with all particles in the right well. 
This is indicated by C in Fig.^. The situation where it is 
possible to have the same values of a chemical potential 
for states with all particles in the left or right well is in 
fact a necessary condition for loop-like structure to be 
present. Recall that in the non-interacting case, match- 
ing of corresponding energy levels is responsible for an 
avoided crossing. 

For 5 > it is possible to have a similar scenario to 
the attractive particles case. The only exception is that 
the ground energy level does not reveal any loop (see 
Fig. 0). Indeed, to transfer the probability density of the 
ground state from the harmonic well to the well created 
by the laser beam we have to move the beam further and 
further towards the center because, contrary to the at- 
tractive particles case, taking the density from the right 
well leads to a deepening of the effective potential (Q) at 
that place. 

Considering the strange shapes of levels in Fig. || one 
may wonder why we have been so successful with time- 
dependent numerical approach. At the beginning of 
the excitation process, the parameter of the system is 
changed sufhciently slowly to follow "adiabatically" a sin- 
gle level far from a loop ("avoided crossing") region. In 
that region, on the other hand, we are diabatic in a sense 
that there is insufhcient time for the wavefunction to 
change its shape appreciably. Thus for non-interacting 
particles we make a diabatic jump from one branch of 
the avoided crossing to the other (following the wave- 
function). Similarly, for interacting particles, loops are 
passed quickly in a way which locally minimizes any sig- 
nificant change of the wavefunction shape. 

To express this more quantitatively we can employ the 
Hellman-Feynman theorem. For the linear Schrodinger 
equation, the slope of the level with respect to a change 
of a parameter, say xq, is simply the expectation value 
(calculated with the help of the corresponding wavefunc- 
tion) of the derivative of the Hamiltonian with respect to 
xq. a diabatic passage does not significantly change the 
wavefunction, so the change of the slope is minimized. 
For the GPE, an analog of the Hellman-Feynman theo- 
rem can be written as 

dfi , ,dV dip dip* 

1— "^ \'p\-i— + 9l—'f +9^-j—m^ (10) 

axf) dxQ dxQ dxo 

and again if the change of the potential is so quick that 
the wavefunction does not react significantly, the slope 
of the level remains the same. 



V. SUMMARY AND CONCLUSIONS 

The proposed method of efficient cohective excitation 
of a BEC seems to be quite robust and aUows preparing 
both sohtonic and vortex-hke excitations of the conden- 
sate. Although the method is complementary to other 
schemes, some of which have been already implemented 
experimentally, still it may be advan tag eous in some 
cases. As we have discussed already in |22| ] our approach 
is somehow closest in spirit to the adiabatic scheme pro- 
posed in ||l^. That approach seems also to be quite ro- 
bust and allows various excitations to be created. The 
method of ||l^ effectively uses two atomic internal states, 
and thus involves a "two-component" condensate. Our 
approach does not entangle internal and external excita- 
tions - this may be advantageous in some applications. 
More importantly our diabatic (or rather "generalized 
diabatic") method takes necessarily a short time - typ- 
ically of the order of a few periods of the trap in our 
runs. This has to be compared with several hundreds of 
periods necessary for the excitation using the adiabatic 
process |16 . 

Needless to say while one may argue about the advan- 
tages of the proposed scheme a most straightforward way 
to test it would be an experimental approach. It seems 
that both the linear potential sweep for cigar shaped BEC 
or the spiral like excitation for disc shaped quasi 2D con- 
densates require relatively minor changes in the already 
existing laboratory set-ups. 

We have found that "dynamics" of chemical potential 
in the time-independent GPE reveals interesting loop- 
like structures. This prevents us from interpreting the 
changes of the chemical potential levels, with respect to 
the parameter, as some form of level dynamics known 
from the linear Schrodinger equation. The observed 
"hysteresis-like" behavior has its origin in the nonlinear- 
ity of the GPE as explained in the text. Its relation to 
the linear quantum mechanics of multiparticle theory is 
being currently investigated. 

To summarize, we have proposed a simple scheme 
which enables us to create collective excitations (both 
solitons and vortices) of the Bose-Einstein condensate. 
This scheme may serve, we believe, as an alternative to 
other proposed methods already utilized experimentally. 

Note added in proof: The scheme considered by us, 
similarly to other works on this subject p5|- [l9| , p4[ , is 
based on the applicability of the GPE to the descrip- 
tion of collectively excited states of the BEC. That has 
been questioned very recently |35[ . 
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